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Atomic clocks based on ultranarrow 3 P — 1 <So tran- 
sition in divalent atoms may offer a new level of time- 
keeping accuracy. In this scheme the atoms are con- 
fined in optical lattice. The lattice laser wavelength is 
selected in such a way that the dominant perturbation 
of the clock frequency, the induced ac Stark shifts, for 
both clock states exactly cancel. Although other effects 
still perturb the clock frequency, estimates [l| indicate 
that the projected fractional uncertainty of such clocks 
may be as low as 10~ 18 . By comparison, 10~ 15 is the 
fractional uncertainty of the current Cs standard realiz- 
ing the SI definition of the unit of time. This apparent 
advantage of the optical-lattice clocks has motivated a 
number of recent proposals: original Katori's scheme |3| 
with fermionic Sr isotopes has been extended to Mg y], 
Ca , and Yb [I| atoms and to bosonic isotopes [y, Q • 
In addition, various schemes of probing the highly- 
forbidden nsnp 3 Pq — ns 2 So clock transition have been 
proposed: three-photon transition, electromagnetically- 
induced transparency, and transition assisted by external 
magnetic field H0|1]. 

Considering advantages of optical lattice clocks, here 
we investigate an important systematic effect of the 
black-body radiation (BBR) on the frequency of the 
3 Pq — 1 So clock transition. Indeed, the SI definition of 
the second explicitly involves atomic clock operating at 
the absolute zero of temperature. In a laboratory envi- 
ronment with an ambient temperature T, one needs to 
introduce the T-dependent BBR correction to the ob- 
served frequency. Here, using techniques of many-body 
relativistic atomic structure, we compute the BBR shift 
for Mg, Ca, Sr, and Yb and evaluate uncertainties of 
the calculations. As summarized in Table [I] the result- 
ing fractional uncertainties in the clock frequencies at 
T = 300 K are large, ranging from 1 x I0 -17 for Mg to 
3 x 10~ 16 for Yb. 

The main conclusions of this paper are (i) the present 
uncertainty in our computed BBR shift is an obstacle on 
the way towards the projected 10~ 18 accuracy goal; (ii) 
due to T 4 scaling of the BBR shift, it may be beneficial 
to operate at low temperatures, e.g., at liquid nitrogen 
temperatures; (iii) if operating at room temperatures, 



high-precision (0.02%-accurate for Sr) measurements of 
the BBR shifts or related quantities are required; (iv) 
Mg-based clock is the least susceptible to BBR; com- 
pared to Sr, the Mg BBR shift is an order of magnitude 
smaller (see Table HJ. Additionally, we develop a general 
relativistic theory of the BBR shift caused by multipolar 
(EJ and MJ) components of the radiation field. 

TABLE I: Black-body radiation shift for clock transitions be- 
tween the lowest-energy 3 Po and ^So states in divalent atoms. 
^bbr is the BBR shift at T = 300 K with our estimated un- 
certainties, uq is the clock transition frequency, and 8vbbr/vo 
is the fractional contribution of the BBR shift. The last col- 
umn lists fractional errors in the absolute transition frequen- 
cies induced by the uncertainties in the BBR shift. 

Atom 5^bbr, Hz i>o, Hz Sisbbh/vo uncertainty 

Mg -0.258(7) 6.55 x 10 14 -3.9 x 10" iB 1 x 10"" 

Ca -1.171(17) 4.54 x 10 14 -2.6 x 10 -15 4 x 10" 17 

Sr -2.354(32) 4.29 x 10 14 -5.5 x 10~ 15 7 x 10" 17 

Yb -1.25(13) 5.18 x 10 14 -2.4 x 10~ 15 3 x 10~ 16 



The paper is organized as follows. Firstly, we derive 
relativistic formulae for the BBR shift which incorporate 
various multipolar contributions of the radiation field. 
Secondly, we apply these general expressions to deter- 
mination of the BBR shifts in optical-lattice clocks with 
Mg, Ca, Sr, and Yb atoms. Unless specified otherwise, 
atomic units and the Gaussian system of electromagnetic 
units are used throughout. In these units, c = 1/cv, where 
a 1/137 is the fine-structure constant. Temperature 
is expressed in units of E^/ks, where Eh is the Hartree 
energy and ks is the Boltzmann constant. 

Multipolar theory of the black-body radiation shift. The 
BBR shift is caused by perturbation of the atomic en- 
ergy levels by the oscillating thermal radiation. Both 
atomic levels involved in the clock transition are per- 
turbed and the overall BBR correction is a difference of 
the BBR shifts for the two levels. We find that deter- 
mining shift for the upper 3 Pq level requires certain care. 
This level is a part of the 3 Pj fine-structure manifold, 
J = 0,1,2. The separation between the levels in the 
manifold is comparable to the characteristic wavenum- 
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ber of the BBR radiation, 208.51cm- 1 , at T = 300 K, 
and contributions of the BBR-induccd magnetic-dipolc 
and electric-quadrupole transitions to the levels of the 
manifold may be enhanced. Taking these induced transi- 
tions into account requires going beyond the conventional 
electric-dipole approximation, Ref. @- 

Considering a potential importance of the multipolar 
contributions, in this section we derive the relevant for- 
mulae for BBR- induced energy shifts. Although we de- 
velop a general relativistic theory including retardation, 
at the end of the section we reduce our expressions to 
practically important non-relativistic non-retarded for- 
mulas. The derived expressions are generalization of the 
electric-dipole BBR shift of Ref. 0. 

The BBR spectral density is given by the Plank for- 
mula 
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It is a weak perturbation and the time-evolution of the 
reference state \g) for off-resonance excitations can be 
computed assuming that the excited state amplitudes 
adiabatically follow that of the reference state. With a 
generalization to narrow resonant contributions |9(, the 
BBR shift is given by 

L J 07? L J dq , L J gp L J pg 
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with averaging over photon propagation directions k and 
polarizations e. The second-order summation is over in- 
termediate atomic states \p) and involves the Coulomb- 
gauge couplings = (a ■ e) exp (k • r)), a encap- 
sulating the conventional Dirac matrices. P.V. denotes 
the Cauchy's principal value; as elucidated in Ref. |9( it 
is required for a proper treatment of resonant contribu- 
tions. 

While evaluating matrix elements of operators ft", 
we use multipolar expansion of e*v k ' r ) in vector spherical 
harmonics and express the resulting couplings in terms 
of traditional multipole moments Q 
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Here A = is for magnetic (MJ) and A = 1 is for electric 
(EJ) multipolar amplitudes. Explicit relativistic expres- 
sions with retardation for matrix elements of Qj^f can 
be found in Ref. [Tof. Neglecting retardation effects (i.e., 
in the long-wavelength approximation) Qj]^ become 



frequency-independent EJ moments Q 
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where Cjm (r) are normalized spherical harmonics. In 
the case of magnetic-dipole transitions in the nonrela- 
tivistic limit = — § (L + 2S) M . Notice that the re- 
tardation brings correction in the order of {aw) 2 to these 
expressions. 

Averaging over polarizations and propagation direc- 
tions in Eq. (JTJ, we find that the BBR shift is a sum over 
multipolar contributions: 5E g = X)/a ^F g J> ^\ 



Here Oig^ (u>) are the generalized dynamic multipolar 
scalar polarizabilitics 
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A cursory examination of these formulas reveals that 
compared to 2 J multipole, the contribution of 2 ,7+1 mul- 
tipole is suppressed by a factor of a 2 rj (1/137) 2 . Also 
for the same J the magnetic contribution is a 2 weaker 
than that of the EJ photons. As in the theory of multi- 
polar radiative transitions E(J+1) and MJ contributions 
are of the same order in a. 

To illuminate the T-dependence of contributions of in- 
dividual intermediate states we recast the BBR shifts into 
a form ( J g is the total angular momentum of the refer- 
ence state, {gWQ^jlj \\p) i s the reduced matrix element) 
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with universal functions 
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Functions Fj (y) are multipolar generalizations of func- 
tion F (y) introduced by Farley and Wing in the El 
case. We plot our computed Fj functions for J = 1, 2, 3 
in Fig. 2] Fj(y) are odd functions with respect to y. 
From examining Fig.^ it is clear that Fj rapidly change 
around y ~ 1 and slowly fall off for y 3> 1. Depending 
on the value of excitation energy, uj pg — yT, a particular 
intermediate state may introduce either negative or pos- 
itive BBR shift. Notice that Fj are broad distributions, 
they have comparable values for \y\ < 20; this will have 
implications for interpreting our results. 

At large values of the argument \y\ 3> 1, Fj(y) oc 1/y. 
The limit y 1 corresponds to the case when the tran- 
sition energy is much larger than T . If all virtual tran- 
sitions satisfy this requirement, then the leading contri- 
bution to the multipolar BBR shift can be expressed in 
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FIG. 1: Universal multipolar functions Fj(y), Eq. JSJ, for 



terms of static polarizabilitics 
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where C is the Riemann zeta-function. As the scaling 
factor, a 2J + 1 T 2J+2 , is expressed in atomic units, we 
observe that as multipolarity J increases by one, in addi- 
tion to the usual a 2 suppression, there is a temperature 
suppression factor of (k B T/E h f. For T — 300 K this 
suppression is sizable, as (fcsT '/ 'Eh) 2 ~ 9.0 x 10~ 7 . 

BBR shift for 3 Po — x So transition in divalent atoms. 
Below we apply the developed formalism to computing 
the BBR shift for the ^o— 3 Po clock transition in divalent 
atoms. We will assume that the atoms are at the ambient 
temperature of T = 300 K. Both clock levels experience 
the BBR shift and the total shift 5^bbr is the difference 
between the two individual shifts, Svbbr = <5^bbr( 3 Po) — 
<5i / bbr( 1 5'o)- 

Consider first the BBR shift of the ground ^So state. 
Here transition energies of various multipolar transi- 
tions to the upper levels are much larger than T, i.e., 
we are in the y 3> 1 limit of Fig. ^ Here compared 
to the dominant El-induced shift, the contribution of 
Ml transitions is suppressed by a 2 ~ 10~ 4 and E2 by 
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Higher-order multipoles are 



suppressed even more. As to the retardation effects in 
El matrix elements, we expect that they would be sup- 
pressed by a factor of a 2 (ksT / Eh) 2 ~ 10~ 10 . Neverthe- 
less, since the fractional contribution of the BBR shift to 
the clock frequency is at 5 x 10~ 15 level (see Table QJ, 
one would need to introduce the Ml corrections at the 
projected accuracy of 10 -18 . 

For the 3 Po levels, the characteristic thermal pho- 
ton frequency is comparable to the fine-structure inter- 
vals for the 3 Pj manifold. The 3 Po level is connected 
by Ml transition to the 3 Pi level and by E2 transi- 
tion to the 3 P2 level. For these transitions the values 
of the relevant functions Fj ~ 1, see Fig. and we 
estimate SE g M1) ~ a 2 (aT) 3 , 8E ( g E2) ~ (aT) 5 , while 



6E, 



■jg E ~ a 3 (T) 4: /LU3 Dl _3p . Our numerical estimate, 
based on the transitions inside the fine-structure man- 
ifold lead to the following values of the BBR shifts for 



Sr: 6E { g M1) = 2.4 x 10~ 5 Hz and 6E {E2) = 2.5 x 10~ 8 Hz. 
Since the El BBR shift for Sr is ~ 2 Hz, the Ml and E2 
contributions can be neglected at the present 1%-level of 
accuracy of our calculations. 

We find that although the thermal photon energy is 
close to the fine-structure intervals, the induced multi- 
pole BBR shifts are not amplified. The main reason is 
that the BBR energy distribution is broad: the universal 
functions Fj have comparable values for a wide range of 
excitation energies, \uj\ < 20 T, see Fig. ^ For example, 
for Sr 3 P - 3 Pi El transition Pi w 0.16, while for the 
3 P - 3 Pi Ml transition F 1 w -0.41 and for the 3 P - 3 P 2 
E2 transition P 2 p» —0.36. For such a broad distribution, 
the multipolar BBR shift is determined by the prefactor 
in Eq. resulting in a suppression of multipoles beyond 
El. 

Based on the above discussion, we may exclusively fo- 
cus on the electric-dipole (J = 1, A = 1) contribution to 
the BBR shift. ^From our general expressions we obtain 
an approximate formula, 
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/T and a g (0) is the traditional static 
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dipole polarizability. To arrive at the above equation, 
we used asymptotic expansion Pi (y) « ||- - 1 - ■^ 2L 
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which has an accuracy better than 0.1% for 
\y\ > 10. 77 represents a "dynamic" fractional correction 
to the total shift. The leading contribution is determined 

( El) 

by polarizability and below we compute a g (0) using 
methods of atomic structure. 

Evaluation of the static dipole polarizabilitics follows 
a procedure of Ref. |ll|. Here we only outline the ma- 
jor steps. We employ relativistic many-body code de- 
scribed in Refs. 0, [r|, 13 • The employed formalism 
is a combination of configuration-interaction method in 
the valence space with many-body perturbation theory 
for core-polarization effects. In this method, one deter- 
mines wave functions from solving the effective many- 
body Schrodinger equation 



{P FC + S(P)} |*„) = E n |*„) , 



(8) 



Here Pfc is the frozen-core Dirac-Hartree-Fock Hamil- 
tonian and self-energy operator £ is a core-polarization 
correction. To improve upon this approximation, one 
can introduce an adjustable energy shift 5 and replace 
E(P) — > H(E — S) in the effective Hamiltonian, Eq. JSJ. 
We have determined S empirically, from a fit of theoreti- 
cal energy levels to experimental spectrum. Inclusion of 
this shift mimics high-order corrections in perturbation 
theory. In addition, we incorporated dressing of the ex- 
ternal electromagnetic field (core shielding) in the frame- 
work of the random-phase approximation (RPA). To find 
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TABLE II: Static electric dipole polarizabilities in a.u. and 
BBR shifts in Hz for the ground ^So and the lowest-energy 
3 Po excited states of Mg, Ca, Sr, and Yb atoms. Theoretical 
uncertainties are indicated in parenthesis. 





Mg 


Ca 


Sr 


Yb 


ai s 0i 

SUBBR.( So) 


71.3(7) 
-0.614(6) 


157.1(1.3) 
-1.353(11) 


197.2(2) 
-1.698(2) 


120.5(3.1) 
-1.04(3) 


asp 

SubBR^Po) 


101.2(3) 
-0.872(3) 


290.3(1.5) 
-2.524(13) 


458.3(3.6) 
-4.052(32) 


266(15) 
-2.29(13) 



valence contribution to a g (0) we summed over the in- 
termediate states in Eq. (0 using the Dalgarno-Lewis- 
Sternheimer method (la ]. A small correction to polariz- 
ability due to core-excited intermediate states in Eq. © 
was computed within the relativistic RPA. 

The results of calculations for the static electric dipole 
polarizabilties for the ns 2 l S$ and nsnp 3 Pq states are 
presented in Table [H] The listed values of the ground- 
state polarizabilities of Mg, Ca, and Sr were obtained 
by us earlier To estimate their uncertainties we 

used the fact that the intermediate state nsnp P± con- 
tributes to the polarizability at the level of 95-97%. For 
calculating the polarizabilities we used the best known 
in the literature values of the (ns 2 ^Sol \D\ \nsnp x Pi) ma- 
trix elements tabulated in 11]. For instance, for Sr 
(5s 21 S ||£>||5s5p 1 Pi)| = 5.249(2) a.u. leading to 0.1% 
error in aig . The uncertainties in the remaining po- 
larizabilities were estimated as a half of the difference 
between two predictions obtained with 5 = and with 
5 determined with the best fit to the experimental en- 
ergies (thus mimicking omitted higher-order many-body 
corrections). The uncertainties in the ground-state po- 
larizabilities range from 0.1% for Sr to 3% for Yb. For 
the 3 Po states the errors range from 0.3% for Mg to 6% 
for Yb. 

With the computed polarizabilities we can find the 
BBR frequency shifts with Eq.Q). The "dynamic" cor- 



rection r) is negligible for the ^So states, but is needed 
for the 3 Pq calculations. Indeed, for the ground state, 
the smallest excitation energy Ei P o — Ei So is equal to 
21698 cm -1 for Sr. At T = 300 K the characteristic 
value of y ~ 100 for all the atoms. By contrast, for 
the 3 Pq clock level, the transitions to the nearby 3 D\ 
level involve smaller energies. For Sr, the relevant en- 
ergy is only 3841 cm -1 corresponding to characteristic 
value of y ~ 20. At this value, the "static polarizabil- 
ity" approximation, F\(y) ~ 47r 3 /(45y), has only a few 
percent accuracy. While evaluating n we find it sufficient 
to truncate the summation over intermediate states at 
the lowest-energy excitation. This "dynamic" correction 
contributes to the BBR shift of the 3 Pq state at 0.1% level 
in Mg, 1% in Ca, 2.7% in Sr, and 0.7% in Yb. Notice 
that since the clock BBR shift is obtained by subtracting 
BBR shifts of the individual levels, the "dynamic" cor- 
rection contributes at an enhanced 5% level in Sr. These 
"dynamic" corrections must be taken into account if the 
BBR shifts are derived from dc Stark shift measurements. 

Finally, we combine the BBR shifts of the individual 
clock levels and arrive at the overall BBR corrections 
summarized in Table [I] Our computed BBR shift for 
Sr, —2.354(32) Hz is in agreement with an estimate pj 
of —2.4(1) Hz. Our uncertainties are better than 3%, 
except for Yb where the error is 10%. As discussed in 
the introduction, although resulting from state-of-the-art 
relativistic atomic-structure calculations, these errors are 
still large and substantially affect the projected 10~ 18 
fractional accuracy of the lattice-based clocks (see Ta- 
ble P). Potential solution involves operating the clocks at 
cryogenic temperatures. 
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